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Abstract The strongly interacting Quark-Gluon-Plasma (sQGP) created in relativis- 
^ , tic nucleus-nucleus collisions, can be described by hydrodynamic models. Low energy 

hadrons are created after the so called freeze-out of this medium, thus their distributions 
reveal information about the final state of the sQGR Photons are created throughout the 
evolution, so their distributions carry information on the whole expansion and cool down. 
We show results from an analytic, 1+3 dimensional perfect relativistic hydrodynamic so- 
lution, and compare hadron and photon observables to RHIC data. We extract an average 
equation of state of the expanding quark matter from this comparison. In the second part 
of this paper, we generalize the before mentioned analytic solution of relativistic perfect 
fluid hydrodynamics to arbitrary temperature-dependent Equation of State. We investigate 
special cases of this class of solutions, in particular, we present hydrodynamical solutions 
with an Equation of State determined from lattice QCD calculations. 

1 Introduction 

The almost perfect fluidity of the experimentally created strongly interacting Quark-Gluon-Plasma 
at the Relativistic Heavy Ion Collider (RHIC) [1] showed that relativistic hydrodynamic models can 
be applied in describing the space-time picture of heavy-ion collisions and infer the relation between 
experimental observables and the initial conditions. 

In this paper we investigate the relativistic, ellipsoidally symmetric model of Ref. [2]. Hadronic 
observables were calculated in Ref. [3], while photonic observables in Ref. [4]. We also show new 
solutions, which can be regarded as generalizations of the model of Ref. [2] to arbitrary, temperature 
dependent speed of sound, originally published in Ref. [5]. 

2 Equations of hydrodynamics 

We denote space-time coordinates by x 11 = (f, r), with r = (r x , r y , r z ) being the spatial three-vector 
and t the time in lab-frame. The metric tensor is g^, - diag{\, -1,-1,-1). Coordinate proper-time 

is defined as r = -\jt 2 - |r| 2 . The fluid four-velocity is u M = y(l, v), with v being the three-velocity, 
and 7=1/ y/l - |v| 2 . An analytic hydrodynamical solution is a functional form for pressure p, energy 
density e, entropy density cr, temperature T, and (if the fluid consists of individual conserved parti- 
cles, or if there is some conserved charge or number) the conserved number density is n. Then basic 
hydrodynamical equations are the continuity and energy-momentum-conservation equations: 

d IJ (nu> l ) = and d Y T^ v = 0. (1) 

The energy-momentum tensor of a perfect fluid is 

T" v = {e + p)uf 1 u v -pff» v . (2) 
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The energy-momentum conservation equation can be then transformed to (by projecting it orthogonal 
and parallel to u?, respectively): 

(e + p) u v d v u» = (</ v - u»u v ) d v p, (3) 
(e + p) d v u v + u v d v s = 0. (4) 

Eq. (3) is the relativistic Euler equation, while Eq. (4) is the relativistic form of the energy conservation 
equation. Note also that Eq. (4) is equivalent to the entropy conservation equation: 

d„ (cru") = 0. (5) 

The Equation of State (EoS) closes the set of equations. We investigate the following EoS: 

s = K(T)p, (6) 

while the speed of sound c s is calculated as c s = ^Jdp/de, i.e. for constant k, the relation c s = 1 / V* 
holds. For the case when there is a conserved n number density, we also use the well-known relation 
for ideal gases: 

p = nT. (7) 

For A-(r) = constant, an ellipsoidally symmetric solution of the hydrodynamical equations is pre- 
sented in Ref. [2]: 

*-£ -4 4 4 • 

where «o and Tq correspond to the proper time when the arbitrarily chosen volume Vo was reached 
(i.e. T = V 1/3 ),andv(s) is an arbitrary function of s. The quantity s has ellipsoidal level surfaces, and 
obeys u v d v s = 0. We call s a scaling variable, and V the effective volume of a characteristic ellipsoid. 
Furthermore, X, Y, and Z are the time (lab-frame time f) dependent principal axes of an expanding 
ellipsoid. They have the explicit time dependence as X = X$t, Y = Ytf, and Z = Zof, with Xq, Yq, Zq 
constants. 



3 Photon and hadron observables for constant EoS 

From the above hydrodynamic solution with a constant EoS, source functions can be written up. For 
bosonic hadrons, it takes the following form [3]: 

S{x,p)cTx = N -; r , (9) 

n(x)exp^p fl u^(x)/T(x)J - 1 

where N = g/(2n) 3 (with g being the degeneracy factor), H{f) is the proper-time probability dis- 
tribution of the freeze-out. It is assumed to be a 5 function or a narrow Gaussian centered at the 
freeze-out proper-time to- Furthermore, [i(x)/T(x) = lnn(x) is the fugacity factor and d 3 Z ll (x)p fl is 
the Cooper-Frye factor (describing the flux of the particles), and d^E^x) is the vector-measure of the 
freeze-out hyper-surface, pseudo-orthogonal to if. Here the source distribution is normalized such as 
J S (x, p)d 4 xd 3 p/E = N, i.e. one gets the total number of particles N (using c-1, H-l units). Note that 
one has to change variables from t to t, and so a Jacobian of dr/dt - t/r has to be taken into account. 
For the source function of photon creation we have [4]: 

, , p u d 3 Z fl (x)dt , . 

S(x,p)d 4 x = N' -— =N' — r d 4 x (10) 

exp (p fl u> 1 (x)/T(x)) - 1 e\p(pf 1 u^(x)/T(x)J - 1 

where p^d 3 ^ is again the Cooper-Frye factor of the emission hyper-surfaces. Similarly to the previous 
case, we assume that the hyper-surfaces are pseudo-orthogonal to u^, thus d 3 E^{x) = u fl d 3 x. This 
yields then p^if which is the energy of the photon in the co-moving system. The photon creation is 
the assumed to happen from an initial time f, until a point sufficiently near the freeze-out. From these 
source functions, observables can be calculated, as detailed in Refs. [3,4] 
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Dataset 




N\ and HBT 


elliptic flow 


N t 






hadrons 


hadrons 


photons 


Central FO temperature 


T [MeV] 


199±3 


204±7 


204 MeV (fixed) 


Eccentricity 


e 


0.80±0.02 


0.34±0.03 


0.34 (fixed) 


Transverse expansion 


uj/b 


-0.84±0.08 


-0.34±0.01 


-0.34 (fixed) 


FO proper-time 


T [fm/c] 


7.7±0.1 




7.7 (fixed) 


Longitudinal expansion 


zlib 


-1.6+0.3 




-1.6 (fixed) 


Equation of State 


K 






7.9 ±0.7 



Table 1. Parameters of the model determined by hadron and photon observables data. See details in Refs. [3,4]. 




Figure 1. Fits to invariant momentum distribution of pions [6] (left), HBT radii [7] (middle) and elliptic flow [8] 
(right). See the obtained parameters in Table 1. 
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Figure 2. Fit to direct photon invariant transverse momentum data [9] (left), comparison to elliptic flow data [10] 
(middle) and direct photon HBT predictions (right). See the model parameters in Table 1. 

4 Comparison to measured hadron and photon distributions 

Observables calculated from the above source functions were compared to data in Refs. [3,4]. Hadron 
fits determined the freeze-out parameters of the model [3]: expansion rates, freeze-out proper-time and 
freeze-out temperature (in the center of the fireball). When describing direct photon data [4], the free 
parameters (besides the ones fixed from hadronic fits) were k (the equation of state parameter) and f,, 
the initial time of the evolution. 

We compared our model to PHENIX 200 GeV Au+Au hadron and photon data from Refs. [6-9]. 
Results are shown in Figs. 1 and 2, while the model parameters are detailed in Table. 1. The EoS result 
from the photon fit is k — 7.9 ± 0.7 stal + l.5 syst , or alternatively, using k — 1 jc\ 

c, = 0.36±0.02 rttf ±0.04 w< (11) 

which is in agreement with lattice QCD calculations [11] and measured hadronic data [12, 13]. This 
represents an average EoS as it may vary with temperature. The maximum value for f, within 95% 
probability is 0.7 fm/c. The initial temperature of the fireball (in its center) is then: 



T; = 507 + I2 stat ± 90™ Jf MeV 



(12) 
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at 0.7 fm/c. This is in accordance with other hydro models as those values are in the 300 - 600 MeV 
interval [9]. Note that the systematic uncertainty comes from the analysis of a possible prefactor, as 
detailed in Ref. [4]. 

Using the previously determined fit parameters, we can calculated the elliptic flow of direct photons 
in Au+Au collisions at RHIC. This was compared to PHENIX data [10], as shown in Fig. 2, and they 
were found not to be incompatible. We also calculated direct photon HBT radii as a prediction, and 
found R out to be significantly larger than R^de- 



5 New solutions for general Equation of State 



We found new solutions to the relativistic hydrodynamical equations for arbitary s — k (T) p Equation 
of State, as detailed in Ref. [5]. These are the first solutions of their kind (i.e. with a non-constant 
EoS). In the case where we do not consider any conserved n density, the solution is given as: 

cr = (To^, (13) 
v? = —, (14) 

T 

^ = exp U (~ + ^)TT^) d ^- (15) 

For the case when the pressure is expressed as p = nT with a conserved density n, another new solution 
can be written up as: 

T 3 

» = «o-r, (16) 



T 3 ' 



— 

T ' 



(17) 

/3d/3 



= ex p(£ ^±[ K{/3}/ 3\]d/3\. (18) 



Quantities denoted by the subscript (n () , T , cr Q ) correspond to the proper-time t , which can be 
chosen arbitrarily. If for example to is taken to be the freeze-out proper-time, then To is the freeze- 
out temperature. These solutions are simple generalizations of the v(s) = 1 case of the solutions of 
Ref. [2], and the latter also represents a relativistic generalization of the solution presented in Ref. [14]. 

It is important to note that the conserved n solution becomes ill-defined, if -^(k(T)T) > is not 
true. In such a case, one can use the solution without conserved n (Eqs. (13)— (15)). 

If k is given as a function of the pressure p and not that of the temperature T, a third new solution 
can be given as: 



if = — , 



- =exp 



OS) + d K (J3) 



48 



48 )k(J3)+1 



(19) 
(20) 

(21) 



i.e. almost the same as in Eq. (15), except that here the integration variable is the pressure p. 
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Figure 3. Left: The temperature dependence of the EoS parameter k from Ref. [1 1] is shown with the solid black 
curve. In the shaded T range (173 MeV - 230 MeV) (k(T)T) (red dashed line) becomes negative, thus the 
solution shown in Eqs. (13)— (15) shall be used with this EoS. Right: Time dependence of the temperature T{t) 
(normalized with the freeze-out time t/ and the freeze-out temperature Tf) is shown. The four thin red lines 
show this dependence in case of constant k values, while the thicker blue lines show results based on the EoS of 
Ref. [11]. 



6 Utilizing a lattice QCD EoS 



A QCD equation of state has been calculated by the Budapest- Wuppertal group in Ref. [11], with 
dynamical quarks, in the continuum limit. In their Eq. (3.1) and Table 2, they give an analytic 
parametrization of the trace anomaly 7 = e- 3/?asa function of temperature. The pressure can also 
be calculated from it, as (if using the normalized values and % = c = 1 units) ^ = y4tt*- F rom 
this, we calculated the EoS parameter k = I/p + 3 as a function of the temperature, as shown in Fig. 3 
(left plot). Since in a T range ^ (k(T) T) becomes negative, the solution without conserved number 
density n (presented in Eqs. ( 1 3)— ( 15)) was used. [5] 

We utilized the obtained k(T) and calculated the time evolution of the temperature of the fireball 
from this solution of relativistic hydrodynamics. The result is shown in Fig. 3 (right plot). Clearly, 
temperature falls off almost as fast as in case of a constant k — 3, an ideal relativistic gas. Hence a 
given freeze-out temperature yields a significantly higher initial temperature than a higher k (i.e. a 
lower speed of sound c s ) would. 

Let us give an example! We shall fix the freeze-out temperature, based on lattice QCD results, to a 
reasonable value of Tf = 170 MeV. Let all the quantities with subscript correspond the the freeze- 
out, we shall thus index them with /. In this case, already at 0.3 x Tf (30% of the freeze-out time), 
temperatures 2.5-3x higher than at the freeze-out can be reached. To give a full quantitative example, 
let as assume the following values: 

Tf = 8 fm/c and T; n j t = 1.5 fm/c, then (22) 
T f = 170 MeV => T m * 550 MeV (23) 

(and even higher if T; n i t is smaller). This value would have been reached with a constant EoS of k « 4, 
even though the extracted average EoS values are usually above this value. The reason for it may be, 
that for the largest temperature range, k obtains values close to 4, as shown in the left plot of Fig. 3. 

In general, the mentioned lattice QCD equation of state of Ref. [11] and our hydro solution yields 
a T(t) dependence. Then, if the freeze-out temperature Tf and the time evolution duration r//r; n i t are 
known, the initial temperature of the fireball can be easily calculated, or even read off the right plot of 
Fig. 3, as it was drawn with units normalized by the freeze-out temperature and proper-time. 
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7 Conclusion 

Exact parametric solutions of perfect hydrodynamics can be utilized in order to describe the matter 
produced in heavy ion collisions at RHIC. We calculated observables from a relativistic, 1+3 dimen- 
sional, ellipsoidally symmetric, exact solution, and compared these to 200 GeV Au+Au PHENIX data. 
Hadronic data are compatible with our model, and freeze-out parameters were extracted from fits to 
these data. [3] 

From fits to direct photon data, we find that thermal radiation is consistent these measurements, 
with an average speed of sound of c s = 0.36 ± 0.02 slat ± 0.04 syst . We can also set a lower bound on the 
initial temperature of the sQGP to 507 ± \2 stat + 90 syst MeV at 0.7 fm/c. We also find that the thermal 
photon elliptic flow from this mode is not incompatible with measurements. We also predicted photon 
HBT radii from this model. [4] 

In the second part of this paper, we have presented the first analytic solutions of the equations of 
relativistic perfect fluid hydrodynamics for general temperature dependent speed of sound (ie. general 
Equation of State). Using our solutions and utilizing a lattice QCD Equation of State, we explored the 
initial state of heavy-ion reactions based on the reconstructed final state. In yJs NN = 200 GeV Au+Au 
collisions, our investigations reveal a very high initial temperature consistent with calculations based 
on the measured spectrum of low momentum direct photons. [5] 
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